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Abstract. We present early results from the ongoing Hydrogen Accretion in LOcal GAlaxieS 
(HALOGAS) Survey, which is being performed with the Westerbork Synthesis Radio Telescope 
(WSRT). The HALOGAS Survey aims to detect and characterize the cold gas accretion process 
in nearby spirals, through sensitive observations of neutral hydrogen (Hi) 21-cm line emission. 
In this contribution, we present an overview of ongoing analyses of several HALOGAS targets. 
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1. Survey description and status 

The scientific motivation for the HALOGAS Survey, and a detailed description of 
the observational setup, are provided by Heald et al. (2011), Briefly, we have selected a 
sample of 24 edge-on and moderately inclined nearby galaxie^ using neutral criteria, 
for deep (120 hr) Hi line observations at WSRT. The observations are sensitive to faint, 
diffuse gas (typical column density sensitivity ~ lO^^cm"^) and to small unresolved Hi 
clouds (typical mass sensitivity of order 10^ Mq). Within the HALOGAS project, each 
galaxy will be carefully studied to ascertain whether signs of gas accretion are present. 
Results from the individual targets will be combined to give the first indication of the 
ubiquity and general characteristics of the cold gas accretion process in nearby spirals. 

At the time of writing, the H I observations are nearly completed; we anticipate that the 
guaranteed WSRT observations (20 of 22 targets) will be finished before summer 2011. 
A number of ancillary programs are also in progress and will supplement the neutral gas 
observations: deep UV (GALEX) and optical (INT) images are being obtained, and a 
novel multi-slit technique will obtain 3D optical spectroscopy of the ionized gas in the 
edge-on survey targets (Wu et al. (2011)). 



t Two of the 24 HALOGAS sample galaxies, NGC 891 and NGC 2403, have already been the 
subject of deep H I observations and are not reobserved in our program. 
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Figure 1. Neutral hydrogen in NGC 1003. Left: Total Hi (moment-0) map overlaid on DSS 
image. Center: Total H i (moment-0) map constructed from HVC analog H i emission only. 
Right: Velocity field (moment-1 map), with cloud positions indicated. 

2. Early results 

First results, based on the first semester of survey observations (the HALOGAS Pilot 



I (2011)1 The survey targets presented in that work 



Survey), are presented by Heald et 

are UGC 2082, NGC 672, NGC 925, and NGC 4565. In this contribution we present 
preliminary analysis for three additional survey targets. A more complete analysis for 
each of these is in preparation. 

2.1. NGC 1003 

NGC 1003 has a star formation rate of OM Mcr, yi-^ (|Heald et al. (2011)|. The width 



of the Hi disk is 21', compared to the roughly 5' optical disk (see Figure [1| . Assuming 
D = llMpc, we obtain a total Hi mass of 5.34 x 10^ Mq. Our tilted-ring modeling 
(described below) indicates a strongly warped, flaring thick disk. Several distinct Hi 
features are detected, as shown in Figure [T] None of the clouds, nor the accretion complex, 
are visible in the optical DSS plates. Cloud 1 has an H I mass of 2 x 10^ Mq. Cloud 2 has 
3.3 X 10^ Mq of H I. Cloud 3 is an unresolved object within the accretion complex with 
an H I mass of 3 x 10^ Mq. The accretion complex, excluding Cloud 3, has an H I mass 
of 2.9 X 10^ M0. The Milky Way and M31 are both known to have intermediate- and 
high-velocity clouds similar to these. The masses we obtain here, as well as the heights 
above the disk, are consistent with those derived by Wakker et al. (2008) [ for Milky Way 
clouds and complexes. Over a dynamical time, they can provide fresh gas compensating 
for only - 2% of the SFR in NGC 1003. 

The GIPSY program GALMOD was used to hand-craft a tilted-ring model of NGC 1003. 
Literature values for scale height, rotational velocity, and other parameters were used, 
and inclination and position angle were estimated. A radial profile was made and used 
for surface brightness values. This model was refined until further improvements could 
no longer be made, and the best parameters were fed into the Tilted Ring Fitting Code 
(TiRiFiC) developed by [Jozsa et al. (2007) | The program fits an arbitrary number of 
parameters at a time for each of the galaxy's rings. We varied a limited number of 
parameters at once (typically two to three) and then fixed them at their optimum values. 
Three final models were produced: a thin-disk model, in which scale height was fixed at 
around 200 pc; a thick-disk model, in which scale height was varied as a unit for all rings 
(and resulted in a scale height of ^ Ikpc), and a flare model where the scale height was 
varied individually for all rings. The east and west halves of the galaxy were modeled 
separately because of an asymmetrical "blob" feature in the eastern side of the disk. 

2.2. NGC 4244 
NGC 4244, having a star formation rate of 0.058 M0yr" 



the low star-forming end of the HALOGAS sample. Its high inclination 



( Heald et al. (2011)[ ), is on 
°) allows 



WSRT HALOGAS Survey 



3 



37°44 




40' 

R.A. (2000.0) 




Right Ascension (J 3000} 



Figure 2. Left: Integrated Hi map of NGC 4244, corrected for primary beam attenuation. 
Contours begin at A^hi = 1-4 x 10^° cm~^ and increase by factors of two. Right: Integrated Hi 
map of NGC 5023 overlaid on a false-color optical picture constructed from DSS images. The 
contours begin at A^hi = 1.9 x 10^^ cm~^ and increase by factors of two. The arrows indicate 
individual extraplanar features of interest. 

for readily determining the vertical H I structure, as well as an accurate assessment of 
the kinematics of the neutral gas. For these reasons it will provide substantial insight 
concerning any connection between star formation and H I halo properties. The total H I 
map is shown in Figure |2j The disk is clearly warped, and does not show evidence for 
a bright diffuse halo component. Through careful tilted-ring modeling (similar to that 
but without the use of TiRiFiC), the disk of NGC 4244 is found to 
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described in 

be warped both parallel to, as well as perpendicular to, the line of sight. The modeling 
confirms that we do not detect a substantially extended H I halo. We do, however, detect 
vertical gradients (lags) in the rotational speed, with magnitudes —12 ± 2kms~^kpc~^ 
and —9 ± 2kms~^kpc~^ in the approaching and receding halves respectively. These 
lags decrease in magnitude to —7 ± 2kms~^ kpc~^ near a radius of 9kpc. The detection 
of the radial variation in the vertical velocity structure will be valuable information for 
comparison with physical models of the disk-halo interface in spiral galaxies. 

Additionally, a few prominent localized features are detected in NGC 4244. Among 
them is a shell in the approaching half, located directly above a region of star formation. 
A second feature in the receding half displays an elongated, curved path in the Hi 
contours, extending away from the major axis to approximately 1.5' above the midplane. 
This feature also appears to correspond with star formation, exhibiting a possible hole or 



vent in the H I, which is partially filled by Ha emission (Hoopes et al. (1999)), indicating 
a possible connection between the two. Finally, a pronged, faint streak is seen over a 
velocity range of ~ 50kms~^ (cf. the rotation speed, 95 — lOOkms"^) in the receding 
half. This feature shows no apparent connection to star formation and exhibits traits 
similar to ram pressure stripping. 

2.3. NGC 5023 

The edge-on galaxy NGC 5023 is a small {v^^i ~ SOkms""*^), slightly warped galaxy at a 
distance of 6.6 Mpc (Heald et al. (2011)| ). Figure [2] shows an integrated Hi (moment-0) 
map. The 21-cm line emission in this galaxy extends on average up to a projected distance 
of ~ 80^' from the plane, with individual features extending to ~ 100'^ (indicated in the 
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figure with arrows). At ttie adopted distance, ttiese vertical distances translate to 2.6 
kpc and 3.2 kpc, respectively. The most vertically extended features also clearly show up 
in the individual channel maps of the data. Even though the bulk of the high-latitude 
emission could be projected above the major axis due to a line-of-sight warp, it is more 
difficult to explain the individual features in this manner. The extraplanar gas also has 
lower projected rotational velocities than the gas in the mid-plane. Detailed modeling 
is in progress to show whether the kinematics is due to a line-of sight warp, a vertical 
gradient in the rotation curve (as seen in e.g. NGC 4244), or some other effect. 



The star formation rate in NGC 5023 is extremely low (0.032 yr~^; Heald et al. 
|(2011)[ ) . Attributing the extended extraplanar H I features to galactic fountain activity 
would therefore seem unlikely. However, deep Ha imaging performed by Rand (1996) | 
shows faint diffuse ionized emission in the extraplanar regions - at the same radial lo- 
cations as the H I features indicated in Figure [2] (they are however not detected to the 
same vertical distances as the Hi features). The correlation between extraplanar diffuse 
ionized gas (EDIG) features and star formation has been well established (e.g. 'RandJ 
|(1996)D . This may indicate that the extraplanar H I features in NGC 5023 are indeed the 
result of star formation activity. Our 3D optical spectroscopic observations (mentioned 
in §[T]) of this target will be very interesting for comparison with the H I data. 



3. Intermediate conclusions and ongoing work 

As the HALOGAS Survey progresses, the number of galaxies with well-studied H I 
morphology and kinematics continues to grow. The extraplanar H I features detected 
in HALOGAS targets have so far been more subtle than, e.g., the massive multiphase 
gaseous halo in NGC 891. However, it is not yet clear which characteristics of the un- 
derlying galaxy dictate the gaseous content of spiral galaxy halos. It does not seem to 
be the case that halos are tied in a straightforward way to the SFR of the host galaxy. 
NGC 4244, for example, has a low SFR and no halo. As shown here, NGC 5023, with 
a similarly low SFR, does have extended extraplanar Hi features. UGC 7321, with an 
even lower SFR, also has an H I halo as shown by Matthews & Wood (2003) 

Which other galactic properties are relevant to the characteristics of gas in halos? 
What role does cold gas accretion play? The full HALOGAS Survey will provide access 
to these questions, by providing uniformly sensitive data over a large sample of galaxies, 
with a broad range in properties such as mass, SFR, and environment. 
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